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In this paper, we study the existence of Floquet topological insulators for VT symmetric non- 
Hermitian Hamiltonians. We consider an array of waveguide in ID with periodically changing 
non-Hermitian potential and predict the existence of Floquet topological insulators in the system. 
We also extend the concept of Floquet topological phase to a two dimensional non-Hermitian system. 


I. INTRODUCTION 


Recent discovery of topological insulators has attracted 
great interest p] (reference therein). A topological insu¬ 
lator has gapped energy spectrum in the bulk while it 
has gapless robust edge states. The concept of topologi¬ 
cal insulator has recently been generalized to periodically 
driven systems HE], Driving a topologically trivial sys¬ 
tem into a non-trivial topological phase is possible by a 
periodic driving. The topological insulators for period¬ 
ically driven systems are known as Floquet topological 
insulator (FTI). Many different groups are searching for 
the existence of both topological and Floquet topologi¬ 
cal insulators in photonics context mm- Haldane and 
Raghu introduced the photonic analogue of topological 
insulator in 2008 [8]. Photonic Floquet topological in¬ 
sulators were first experimentally realized by assuming 
that the periodicity in time is replaced by a periodicity 
in propagating direction [9 • We emphasize that there are 
some differences between the description of standard and 
Floquet topological insulators. For example, the Chern 
number, which provides the net number of chiral edge 
states crossing the band gap, is not good enough to give 
a full characterization of the topological properties of pe¬ 
riodically driven systems. It was shown in m that so- 
called anomalous edge phenomena appears only for FTI. 
The existence of topological phases has recently been 
investigated for non-Hermitian Hamiltonians jT9H25] . 
Of special importance is VT symmetric non-Hermitian 
Hamiltonians, where V and T operators are parity and 
time reversal operators, respectively. This is because of 
the fact that Hermiticity requirement was shown to be 
replaced by the analogous condition of VT symmetry for 
the reality of the spectrum [26] . It is well known that a 
VT symmetric non-Hermitian Hamiltonian admits real 
spectrum as long as non-Hermitian degree is below than 
a critical number. The search for topological phase was 
first considered by Hu and Hughes m and Esaki et al. 
[20]. They found that topological states are not stable 
because of the existence of complex energy eigenvalues. 
In [T9], they considered Dirac-type non-Hermitian Hamil¬ 
tonians and concluded that the appearance of the com¬ 
plex eigenvalues is an indication of the non-existence of 
the topological insulator phase in non-Hermitian models. 
Similarly, topological phase was shown to be unstable for 
a non-Hermitian generalizations of the Luttinger Hamil¬ 
tonian and Kane-Mele model [20], a one dimensional non- 


Hermitian tight binding model [22] and non-Hermitian 
Su-Schrieffer-Heeger (SSH) model with two conjugated 
imaginary potential located at the edges of the system 
m- Recently, we have found stable topological phase 
for a non-Hermitian VT symmetric system for the first 
time in the literature [25] . 

The existence of topological phase transition for VT sym¬ 
metric non-Hermitian systems is an important question. 
In the literature, FTI for a non-Hermitian system has 
not been studied yet. In this paper, we study the exis¬ 
tence of FTI for some non-Hermitian Hamiltonians. We 
first consider ID SSH model with periodically changing 
non-Hermitian potential and show the existence of FTI. 
Secondly, we study a 2D system and predict FTI in the 
2D system. 


II. FLOQUET TOPOLOGICAL INSULATOR 

Consider a non-Hermitian periodical Hamiltonian 
in the propagation direction with the period T, 
H(z ) = H(z T). The energy spectrum of a period¬ 
ical Hamiltonian can be found using the Floquet the- 
ory |2Z|. The Floquet spectrum is 2i t/T periodic in 
quasi energy just as a repeated zone scheme in quasi 
momentum for conventional band structure. Accord¬ 
ing to the Floquet theorem, there exist solutions to the 
Schrodinger equation of the form = e~ l€at & a , where 
is the periodical Floquet eigenstate with period T 
and the time-independent eigenvalues e a are the Floquet 
quasi energies. Note that the quasi energies are con¬ 
stants in time. Let us define the Floquet Hamiltonian 
as Hp = H — id/d z that satisfies Hp<$> a = e& a . Since 
the Floquet states are periodical, they satisfy the rela¬ 
tion <fr an = & a e inujz , where the orthonormality condi¬ 
tion reads << $cm|$/3m >>= ^ Jq dz < ® a n\®/3m > = 
dapfimn and the second bra-ket notation here is used 
to denote the integration over the propagation direction 
z. Of particular interest for a non-Hermitian periodical 
Hamiltonian is the situation when all quasi energies are 
real. In the present paper, we numerically find quasi en¬ 
ergy spectrum for our non-Hermitian Hamiltonian. 
Consider now a ID tight binding waveguide array. Let 
the operators T n and a n create and annihilate a particle 
on site n, respectively. We are interested in the system 
with staggered coupling coefficients (i.e. alternate be¬ 
tween weak and strong coupling coefficient). Let t\ and 


2 


t 2 denote alternating coupling coefficients. This model, 
known as SSH model [28], is a two band model exhibiting 
nontrivial topological properties. In addition to the SSH 
model, suppose periodically changing non-Hermitian im¬ 
purities are inserted in the whole lattice. Let the parame¬ 
ter 7 denotes the non-Hermitian degree of the impurities. 
The impurities with 7 > 0 and 7 < 0 represent gain and 
loss materials, respectively. The Hamiltonian for our sys¬ 
tem reads 

H = — Y t n a^a n+ 1 + h.c. + (-l) n i 7 cos (uz)al l a n (1) 

n 

where the position dependent coupling coefficients t n = 
£1 when n is odd and t n = £ 2 when n is even and uj is the 
angular frequency of the modulation. This Hamiltonian 
describes the propagation of an optical field in a tight 
binding waveguide array with periodically changing gain 
and loss in the propagation direction. If we replace the 
propagation direction z by time £, then the same Hamilto¬ 
nian describes a tight binding bichromatic optical lattice 
with time dependent gain/loss terms. 

We can rewrite the Hamiltonian as H = aJ l iL nm aJ n , 

n,m 

where H nm are the corresponding matrix elements. Pro¬ 
vided the system has translational symmetry (the system 
is subjected to periodic boundary conditions), the matrix 
form of the Hamiltonian can then be written in momen¬ 
tum space as 

H(k) = £2 sm(2ka)<j y — (£1 + £2 cos (2 ka)) a x -\-ij cos(ujz)<j z 

( 2 ) 

where 0 are Pauli matrices, the crystal momentum k runs 
over the first Brillouin zone, —tt/ a < k < tt/ a and a is the 
lattice spacing. Below we will show that the system sup¬ 
ports topologically nontrivial edge states if u 7 ^ 0. Let us 
define t\ = 1 — A cos <I> and £2 = 1 + A cos <F, where A < 1 
is a real valued constant and the phase parameter <f> is an 
additional degree of freedom. To study topological phase 
in our system, we vary <F from 0 to 27r continuously and 
see if the band gap closes during the continuous change. 
Let us first analyze the case with u = 0. In this case, 
the corresponding band structure can be readily found. 

It is given by E = Ty^i + £2 + 2 £i £2 cos (2 ka) — 7 2 . As 
can be seen, the system has two bands. Consider first 
Hermitian limit, 7 = 0 . In this case, the band structure 
is gapped when $ = 0. The band gap decreases by vary¬ 
ing <f> and vanishes at <F = 7t/2. If we further increase 
<f>, the band gap opens again. As a result, topological 
phase transition occurs at <f> = tt/2 (also at = 37t/2). 
Consider now the presence of non-Hermitian term in the 
Hamiltonian, 7 7 ^ 0. In this case, topological states are 
not stable since complex energy eigenvalues appear when 
the phase transition occurs. This can easily be seen 
by setting t\ = £2 and ka = tt/2 on the above energy 
expression. As a result, we say that no stable topo¬ 
logical phase exists in our non-Hermitian system when 
uj = 0. Let us now search for the existence of Floquet 
topological phase by considering periodically modulated 


non-Hermitian impurities, uj 7 ^ 0. In this case, we apply 
Floquet theory to find the quasi energy spectrum. In 
practice, a truncated Floquet Hamiltonian is used in our 
numerical computation. The corresponding matrix is a 
2(2N F + l) x2(2N f + 1) matrix, where N F is large enough 
that the result doesn’t depend on N F [29] . Strikingly, we 
numerically find stable topological phase provided oscil¬ 
lation frequency uj is large and the non-Hermitian degree 
7 is below than a critical value. In the Fig.l, we plot the 
quasi energy spectrum as a function of ka when uj = 2tt 
and 7 = 1 . Note that the spectrum is cj-periodic in quasi 
energy and 27r/a-periodic in crystal momentum k. Just 
as in the Hermitian case, topological phase transitions 
can take place by varying the coupling coefficients as long 
as generic symmetries of the Hamiltonian are preserved. 
Suppose £2 is adiabatically changed from 0.5£i to 1.5£i, 
where t\ = 1 . As can be seen from the figure, the con¬ 
duction and valence bands are separated by a gap when 
£2 < £i- Quantum phase transition occurs at £2 = £1 
since the band gap closes at ka = ^ 7 t/ 2 . If we further 
increase £2 such that £2 > £1 the band gap opens again. 
Two periodically driven systems are defined to be topo¬ 
logically equivalent if the quasi-energy band structures 
can be smoothly deformed from one to the other without 
closing the band gap. So, we say that the two systems 
with £2 < £1 and £2 > £1 are topologically distinct since 
the quasi-energy gap doesn’t stay open during the con¬ 
tinuous deformation that respects the symmetries of the 
Hamiltonian. Our non-Hermitian system is topologically 
trivial when t\ > £2 and nontrivial when t\ < £ 2 . We 
would like to emphasize that the energy spectrum is al¬ 
ways real during the continuous deformation. So the sys¬ 
tem is topologically stable as opposed to the non driven 
system with uj = 0. The topological distinction between 
the two cases t\ < £2 and t\ > £2 leads to the pres¬ 
ence of topologically protected edge states for the for¬ 
mer case. More precisely, both cases are insulators in 
the bulk, but the case with t\ < £2 possess topologically 
protected edge states. To see the edge states, the non¬ 
trivial topological insulator can be placed in a vacuum, 
which is topologically equivalent to a conventional insula¬ 
tor. Let us now study topological edge states for a system 
with N sites. We adopt open boundary conditions with 
n = 1 and n = N being the two edge sites. Note that 
the translational invariance of the system is broken for 
open boundary conditions. We numerically find the cor¬ 
responding quasi energy spectrum. The last figure in the 
Fig.l plots the quasi energy spectrum as a function of <f> 
when N = 40, 7 = 1, A = 0.4 and uj = 27r. The system 
is topologically trivial (nontrivial) when 7 t/2 < 4> < 3tt/2 
(otherwise) since t\ > £2 (£1 < £ 2 ) in this region. In the 
absence of non-Hermitian potential, doubly degenerate 
zero energy states are the topological edge states. The 
presence of non-Hermitian term lifts the degeneracy and 
two symmetrical states appear on both neighboring sides 
of zero energy line. Let us now study evolution of the 
edge states. We start with an initial edge state obtained 
in the limit 7 = 0 . Then we numerically find the evolu- 
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FIG. 2: The propagation of an initial edge state when N = 40. 
The edge state is protected in the topologically nontrivial case 
(a), while it is delocalized in the topologically trivial case (b). 
The density would grow without bound when uj — 0 while 
oscillates if the quasi energy spectrum is real when uj ^ 0 as 
can be seen in (a). 



FIG. 1: The first three figures plot quasi energy spectra of the 
translationally invariant system for three different coupling 
coefficient t^. The last figure plots the quasi energy spectrum 
for the open boundary condition as a function of <f>, where 
we define for t\ — 1 — 0.4cos(<f>), t 2 = 1 + 0.4cos(4>). The 
parameters are given by N — 40, uj = 2tt and 7 = 1. 


tion by assuming that 7 = 1 and uj = 27r. The Fig.2.a 
and Fig.2.b are for topologically nontrivial system with 

= 0 (ti = 0.4, £2 = 1.6) and topologically trivial sys¬ 
tem with = 7 t /2 (tg = 1.4, = 0.4), respectively. The 

edge state is protected in the former case while it is de¬ 
localized in the latter case as expected. As opposed to 
topological edge states for a Hermitian Hamiltonian, the 
density of the edge state in our non-Hermitian system os¬ 
cillates while propagating as can be seen in the Fig.2.a. 
We note that oscillation of the amplitude of the edge 
state is the direct result of the quasi energy spectrum of 
the non-Hermitian Hamiltonian. To this end, we would 
like to emphasize that we numerically see that complex 
energy eigenvalues start to appear if the non-Hermitian 
degree 7 exceed a critical value at fixed uj = 2tt (the criti¬ 
cal non-Hermitian degree 7 c = 3.6 when uj = 2i r). In this 
case, the system is topologically unstable. The real part 
of the spectrum changes dramatically once the complex 
energy appear in the system. 

We have studied the FTI in ID. Now, our aim is to 
study the FTI in a 2D system. Consider a 2-dimensional 
translationally invariant and periodically driven system 
with the period T = 27t /uj 

H = cos(ujz)a z + A (sin (k x a) a x + sin (k y a)a y ) 

+ (A + cos(k x a) + cos(k y a)) cr z (3) 

where A is a real valued constant defined in the interval 
[— 2 , 2 ] and we set A = 1. The system is supposed to 
be translationally invariant and wave numbers k x and k y 
are therefore good quantum numbers. The corresponding 
energy eigenvalues can be found analytically if the system 


is not periodically modulated, uj = 0 , 

E = =F(sin 2 k x x-\- sin 2 k y y+(A+cos k x x+cos k y y) 2 +5E) 1 / 2 

(4) 

where SE = — 7 2 + 2 i 7 (A + cos k x x+cos k y y). In the Her¬ 
mitian limit, 7 = 0 , the gap in the spectrum can be closed 
and reopened by tuning the parameter A. If A = — 2 
and A = 2 the energy gap closes when k x = k y = 0 
and k x = k y = =F?r, respectively. If A = 0 , the energy 
gap closes when k x = 0 , k y = =F?r and k x = =pL k y = 0 . 
For all other values of A the spectrum is gapped, and 
thus nontrivial topological phase exists in the system. In 
the presence of non-Hermitian potential, 7 7 ^ 0, the sys¬ 
tem becomes unstable since the system has complex en¬ 
ergy eigenvalues. Therefore, no stable topological phase 
is available for the non-driven case just as for the ID case 
considered above. Let us now study the quasi energy 
spectrum for the periodical system uj 7 ^ 0. As in the ID 
case, the quasi energy spectrum is real if uj is large enough 
for a given 7. In other words, no complex quasi energy 
eigenvalues appear during the quantum phase transition. 
The Fig.3 plots numerically calculated quasi energy spec¬ 
trum as functions of k x and k y for three different values 
of A at uj = 6tt and 7 = 1 . One can see the band gap 
closing and reopening when we increase A from —0.5 to 
0.5. The system is an insulator when A = =f 0.5 while the 
band gap closes at A = 0. So, we conclude that the two 
cases with A = tO- 5 are topologically distinct. Since 
the quasi energy spectrum is real and quantum phase 
transition occurs during the continuous deformation of 
the Hamiltonian by varying the parameter A, the sys¬ 
tem is said to be topologically stable. Let us now study 
edge states. The vacuum is topologically trivial so plac¬ 
ing the system in the vacuum leads to the edge states. 
Suppose the lattice is infinite along x while there are N 
sites in y axis. The system is subjected to open boundary 
condition along y. Translation invariance holds along x 
and k x is still a good quantum number. The last figure 
in the Fig.3 plots the corresponding quasi energy spec¬ 
trum as a function of k x when A = —0.5, uj = 67 r and 
7 = 1. As can be seen, the spectrum is gapped with gap¬ 
less edge states crossing between valance and conduction 
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FIG. 3: The quasi energy spectra for different values of A. 
The gap closes and reopens by tuning A. The last figure plots 
the quasi energy spectrum for the open boundary condition 
as a function of k x . The parameters are given by N m 20 unit 
cells wide in the y direction and infinite along the x direction, 
uj — 67 r and 7 = 1 . 


bands. The gapless edge states leave their bulk bands 
from nearly k x = T^/4 and cross near k x = 0. These 
states are edge states since their wave functions are con¬ 
centrated at the edges. However, the states in the bulk 
band are extended states. One can notice that the pic¬ 
tures for Hermitian and non-Hermitian cases are similar 
as long as the spectrum of the non-Hermitian Hamilto¬ 
nian is real. The real part of the spectrum changes dras¬ 
tically if the imaginary part of it is different from zero. 
Similar to the ID case, the quasi energy spectrum be¬ 
comes complex if 7 exceed a uo dependent critical value. 
To sum up, we have studied V and T symmetric non- 
Hermitian Hamiltonians in ID and 2D. We consider two 
band model and show that quantum phase transition oc¬ 
curs in our system. We discuss that topological states are 
not stable if the system is not periodically modulated. By 
periodically modulating the non-Hermitian potential, we 
show that it is possible to get stable Floquet topological 
phase. 
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